Generation of graph-state streams 
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We propose a protocol to generate a stream of mobile qubits in a graph state through a single stationary parent 
qubit and discuss two types of its physical implementation, namely, the generation of photonic graph states 
through an atom-like qubit and those of flying atoms through a cavity-mode photonic qubit. The generated 
graph states fall into an important class that can hugely reduce the resource requirement of fault-tolerant linear 
optics quantum computation, which was previously known to be far from realistic. In regard to the flying atoms, 
we also propose a heralded generation scheme, which allows for high-fidelity graph states even under the photon 
loss. 

PACS numbers: 03.67.Bg, 03.67.-a, 42.50.Dv, 42.50.Ex 



The graph state, marked by having peculiar correlations UJ 
and being a universal resource for quantum computation H, 
has been recognized as one of the most important classes of 
many-body entangled states. In the last decade, numerous 
studies have been conducted to generate graph states in var- 
ious systems fs'-^Ts'l and the proof-of-principle experiments 
have been performed with realizing ineasurement-based logic 
operations using graph states lllOl4l3ll . 

An important feature of the graph state is that it can be 
represented by a graph with each vertex representing a qubit. 
For a given graph (specifically, a simple graph 1 14]), the cor- 
responding graph state is defined as the state generated by 
preparing every qubit in state |+) = + |1)), where |0) 

and |1) are the computational basis states, and subsequently 
performing a CPHASE operation between every pair of qubits 
connected by an edge. This definition is very suggestive in 
that one can essentially build up graph states in a systematic 
way, i.e., by adding qubits one by one, connecting two graph 
states, and so forth. This building-up procedure can be per- 
formed even with nondeterministic CPHASE gates 111511 . leav- 
ing us a problem of finding an optimal strategy lfT6ll . 

In most physical situations, however, such graph-state gen- 
erations are largely restricted both spatially and temporally. In 
particular, when mobile qubits, e.g., photons or flying atoms, 
are of major concern, such restrictions are more crucial unless 
massive reroutings and clever spatial allocations are permitted 
even though the mobile-qubit graph states can be useful for 
various quantum information and communication protocols. 
It would thus be crucial towards useful quantum information 
processing to identify the capabilities of individual systems as 
a source of graph states. 

In this paper, we propose a protocol to generate graph states 
of mobile qubits in a highly restricted but commonly faced 
situation wherein a linear stream of mobile qubits come out 
sequentially from a place containing a single stationary qubit, 
which we will call a parent qubit. In particular, we consider 
two types of systems: firstly, a stream of single photons gener- 
ated recuiTently by an atomic (or atom-like tlTIl ) qubit trapped 
in a cavity, and secondly, a stream of flying atoms passin g se - 
quentially through a cavity containing a single photon lUSll . 



For convenience, let us call the former System I and the latter 
System II. 

Experimentally, only pairwise entanglement generation in 
System I and II using a parent qubit has been realized lil8i[T9tl . 
While photonic graph-state generation in System I has been 
considered theoretically [3, 8], the same for the atoms in Sys- 
tem II has not been thoroughly studied. Our protocol offers 
a versatile and structured way of generating a rich variety of 
graph states and further envisions the possibility of the genera- 
tion in System II. It turns out that the graph states correspond- 
ing to a linear chain of arbitrary star graphs iH . as shown in 
FIG. ^a), can be generated. The graph states and variations 
considered in Ref. for System I fall into the same class 
up to local unitary transformations, albeit not explicit therein. 
This type of graph states is indeed of significant importance 
in fault-tolerant linear optics quantum computation (LOQC) 
||20[ I21I1 . Although one can obtain fault-tolerance thresholds 
for LOQC, its requirement of resources — the number of two- 
photon entangled pairs, time steps, and parallel operations — is 
unrealistically high as it heavily relies on off-line preparation 
of particular graph states that is done with an extremely low 
success rate, spending a huge amount of the resources. Those 
graph states that the off-line preparation is aimed at are in fact 
of the same shape as that in FIG. [TJa)- Our scheme would 
thus hugely reduce the resource requirement of fault-tolerant 
LOQC. In regard to System II, experimental errors may be 
due to the loss of the photon that should be kept in the cavity 
throughout the whole process as a parent qubit. In order to 
circumvent this problem, we propose a heralding scheme that 
can single out successful events when the atoms are in high- 
fidelity graph states. The schemes considered in this paper are 
within cuiTent technology. 

Our protocol makes use of two unit operations to gener- 
ate graph states. The first one is what we call a branch- 
ing operation, which adds a vertex stemming from the par- 
ent qubit, as in FIG. [Ub). This operation brings the parent 
qubit p into a combination of the parent and a mobile qubit 
and maps the state of the parent qubit into the combined state 
as follows: |0)p |0)p|+)i and |l)p \l)p\—)i, where 
|— ) = -^(|0) — |1)) and the subscript i denotes the i-th 
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(a) Linear chain of star graphs 
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(b) Branching 
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(c) pulling-out 



FIG. 1: (a) An example of a graph state that can be generated by 
combining the two operations in (b) and (c). (b, c) Two unit oper- 
ations for adding the i-th qubit. The shaded vertex represents the 
parent qubit. 



mobile qubit. This resembles a CPHASE operation between 
the parent and i-th mobile qubits prepared in the |+) state. 
The other operation is what we call a pulling-out operation, 
which adds a vertex taken by the parent qubit while a new 
mobile qubit takes the previous place of the parent qubit, as in 
FIG. [lie). This operation maps the state as |0)p — > |0)i|+)p 
and \l)p — > )p, which is identical to a branching op- 

eration followed by a SWAP operation. While the general 
SWAP operation is a nonlocal operation, it is important to 
note that this particular SWAP operation can be replaced by 
local Hadamard operations. It is easily seen that by combining 
branching and pulling-out operations, one can generate a state 
corresponding to any linear chain of star graphs, as shown in 
FICEa). 

In System I where photons are mobile qubits, single-qubit 
operations are tractable for both the parent and the mobile 
qubits. The essential part is thus the branching operation, as 
local Hadamard operations transform the branching operation 
into the pulling-out operation. For the branching operation, 
one needs to generate a single photon in such a way that an ar- 
bitrary superposition state of the parent qubit a |0)p + /3 
is transformed to a \0)^ + (3 \1) ^ lo"-), where \a±) are 
two orthogonal polarization states of a photon. This can be 
done by using (or slightly modifying) well-studied existing 
schemes for single-photon generation |22]. For example, one 
can realize the branching operation using two ground hyper- 



fine levels |0) 



= -5) and |1> 



i) with total 



spin F = i, as shown in FIG. |2a), to represent a parent 
qubit. As the transitions with Am = ±1 are coupled to the 
cr±-polarized modes of a cavity, Raman transitions can occur 
for both the ground levels by applying a vr-polarized classical 
field, which generates a single cavity photon while flipping 
the atomic qubit. The single photon leaking out of the cavity 
then has a polarization depending on the atomic initial state, 
as we desired. The spin flip of the atomic state is easily cor- 
rectable by a single-qubit operation, or this correction can be 
simply put off to a later time by employing the Pauli frame 
ll23ll . The spontaneous emission can be suppressed by having 
a sufficiently large detuning. As for the time scale, using this 
conventional method of adiabatic transfer (see Ref. ll22ll ). the 
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FIG. 2: Involved atomic levels and transitions for System I (a) and 
for System II (b). Two ground hyperfine levels represent a qubit. 



width of a single-photon pulse can be estimated as some large 
constant (say, ^10) divided by the cavity decay rate, hence 
the repetition rate can be typically as high as tens of kHz. 

Differently from System I, realizing the branching and 
pulling-out operations is not straightforward in System II as 
it is hard to perform single-qubit operations on the parent 
qubit, which is encoded in a cavity-mode photon. In this 
case, the CPHASE and SWAP operations should be directly 
implemented. For this, we consider a A-type level structure as 
shown in FIG.|2|b). The two transitions are coupled, respec- 
tively, to the two orthogonally polarized modes of the cavity, 
but are largely detuned from the resonance. Consequently, 
when the atom just passes through, no atom-cavity interaction 
takes place. In order to turn on the interaction, we control the 
two detunings separately by means of an ac Stark shift. As 
will be shown later, when only one transition is brought into 
resonance, a CPHASE operation can be performed, and when 
both transitions are brought into resonance, a SWAP opera- 
tion can be performed. We can thus perform these operations 
selectively while the flying atom is passing through the in- 
teraction region of the cavity by applying timely pulses that 
induce the ac Stark shifts. Note that as the flying atoms en- 
ter the cavity with a random timing in experiments, the ex- 
istence of an atom should be detected optically for a correct 
timing of the pulses. Although this protocol would benefit 
from the longer coherence time achievable in state-of-the-art 
microwave cavity QED technologies |1.8f], there are still prac- 
tical shortcomings of them that make our discussion for Sys- 
tem II more relevant to optical cavities |24, 25] at the present 
time. For example, while atoms with the required level struc- 
ture are readily available in the optical regime, this is not the 
case in the microwave regime (i.e., Rydberg atoms). 

Whereas the essential component of System I — coherent 
generation of a single photon by a Raman transition — has 
been well established [22], System II needs further analysis, 
which we focus on in the remainder of this paper. As shown in 
FIG. |2|b), the atom has two degenerate ground levels |^) and 
I r) , which represent a qubit, and an excited level | e) . The tran- 
sition 1^) o |e) (|r) o |e)) is coupled to the L(i?)-polarized 
mode of the cavity with coupling strength (gB.) and de- 
tuning Ai (Aji). In case A^ = A^^ = 0, the interaction 
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Hamiltonian can be written as 



H=L,B. 



(1) 



where denotes the annihilation operator for the /i-polarized 
mode of the cavity, (t| = |e) (/|, and a\f = \e) {r\. We assume 
9l = dR = g when the transitions are resonant. If the atomic 
transition is largely detuned (A^ ^ g^), the effective cou- 
pling between the atom and the external field can be regarded 
as being zero, i.e., — gn = 0. 

We first outline our scheme in an idealized situation where 
both the atomic and the cavity decays are absent. At the begin- 
ning, the cavity should be loaded with a single photon. This 
can be done in combination with an entangling operation. For 
this, the first atom is initially prepared in state |e) and a classi- 
cal pulse that induces the ac Stark shift is applied to adjust the 
detunings to achieve — A^ — 0. After straightforward 
calculation based on the Hamiltonian Hi we can find that if 
the pulse duration tq is chosen to be gro — 7r/2\/2, the state 
is transformed to -^(10 + l-R)), which is the two-qubit 
graph state. 

From the second atom, either a branching or a pulling-out 
operation is performed as outlined before. In order to perform 
a CPHASE operation, only Aji is adjusted to be zero while 
Al ^ gL during a period of time ti such that gRTi = tt. This 
operation transforms the state as \r) \R) — > —\r)\R) while 
leaving |/) \L), \l) \R), and \r) \L) unchanged, which is iden- 
tical to a CPHASE operation. In order to perform a SWAP op- 
eration, we adjust both A^ and A/? to be zero during a period 
of time T2 such that gT2 = tt/V^- Solving the Schrodinger 
equation for the Hamiltonian it can be shown that this op- 
eration transforms the state as follows: \l) \L) —\r)\R), 
10 \R) ^ 10 \R), \r) \L) ^ \r) and \r) \R) ^ \l) \L). 
This is identical to a SWAP operation followed by a ZX op- 
eration acted on both qubits, where Z and X are the Pauli 
operators. The correction of these Pauli operations can be de- 
ferred by employing the Pauli frame ll23ll . 

In real experiments, the above idealized situation is not pos- 
sible to achieve due to the finite finesse of the resonator, to- 
gether with the finite lifetime for the excited state of the three- 
level atoms. To model the evolution of the chain of atoms 
as they cross the cavity, we can use a Lindblad equation un- 
der the Markovian approximation for the density matrix of the 
system i26ll 

d Q 

~ X! i^i^t^S + S^l^f' ~ 2(T^£»(t1;) . (2) 

fj.=L,R 

Here and 7^ stand for the cavity damping and the atomic 
spontaneous emission rates for the /x— circularly polarized 
mode, respectively. In writing the master equation (|2]i we have 
assumed that the cavity field is coupled to a bosonic environ- 
ment at zero temperature. 

Clearly, the performance of our protocol depends on the 
ability of the resonator to keep the excitation for a time long 
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FIG. 3: (a) Fidelity and (b) probability of the heralded generation 
of a 3-qubit cluster state versus the cavity damping and spontaneous 
emission rates, (c) Fidelity and (d) probability of the heralded gen- 
eration of the same state versus the cavity damping (here 7^ = kr) 
in the protocol described (circles), and assuming that there is up to 
a 10% uncertainty (uniformly sampled) in the selection of the in- 
teraction time (squares). We take = kr and 7^ = 7^ in all 
cases. Cavity damping and spontaneous emission rates are measured 
in units of the coupling constant, g, which is set equal to unity. 



enough to perform all the operations needed. This means that 
the fidelity of the state generated will be very much affected 
by the damping and spontaneous emission rates. In spite of the 
photon loss, however, the fidelity can be drastically improved 
nearly as high as unity by adding to the original scheme a 
heralding process based upon the detection of the photon leak- 
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ing out of the cavity. Note that the number of excitations 
S/j=L R^ji^iJ' + (^1 ^'^ ideal case is always kept 

to be one, while the environmental effect only decreases it. 
Consequently, once a photon is detected, it is guaranteed that 
all the preceded operations have been performed along with a 
single excitation kept in the system, and hence the fidelity of 
the final state should be high. We will take the state only when 
the photodetector clicks (in the computational basis) after the 
complete set of gates have been applied to the chain of atoms. 

Another important advantage of employing the heralding 
process is that we can deduce the fidelity and the success prob- 
ability for an arbitrary number of atoms in a pseudo-analytic 
way, avoiding to calculate the whole evolution for all the el- 
ements of the huge density matrix. For convenience, let us 
change the notation to account explicitly for the number of 
photons in the cavity for each polarization, i.e. |10) = \L), 
1 01) = \R), while |00) will denote the absence of the pho- 
ton in any polarization. To calculate the evolution of one 
atom together with the cavity field (within the one or zero- 
excitation subspace), we can expand the density matrix us- 
ing the basis {\110) , |;01> , |rlO) , |r01) , |eOO) , \m) , \rOQ)} 
to solve the master equation (|2|i. After the ideal evolution for 
a CPHASE and a SWAP, only the 4 x 4 sub-matrix at the top 
left corner for those states containing one photon will have 
non-zero elements. If the heralding process is assumed, only 
these elements indeed make contributions to the final state. 
It turns out that even with the photon loss, this sub-matrix, 
when renormalized, remains as a pure state, which can be 
written as a |10) + /? |01) |(/'_r), where is a normal- 
ized state representing other portion of the state. For the ensu- 
ing step, this pure state can be taken as the initial state, where 
the evolution can be calculated in the same manner as if a 
state a |10) + /3 |01) was the initial state. This is because the 
interaction exists only between the cavity photon and the new 
atom, while all the earlier atoms left the cavity do not take part 
in the new evolution. 

As a particular example, we have performed numerical 
calculations for the generation of a 3-qubit graph state. In 
the ideal case this linear graph state will be equiva- 



lent to a 3-qubit GHZ state. Fig. [3la) shows the fidelity 
of the state obtained after heralding and the probability for 
this to happen with respect to the cavity damping rate and 
the atomic spontaneous emission rate. As we are comparing 
pure states, |(/)3) and the fidelity can be calculated as 

T = \{<f>3\(f'3'')\^ ■ Typical values of ((?,7, K)/27r are around 
(16,3,1.25) MHz for state-of-the-art experiments with Rb 
atoms [27] or (34,2.6,4.1) MHz for Cs [24]. Taking these 
parameters, the fidelities are T ~ 0.886 and w 0.998, re- 
spectively, while the probabilities of heralding sitV ~ 0.84% 
and V ~ 0.19%. This is relevant in order to calculate the rep- 
etition rate which can be obtained with this protocol. In the 
evolution we have also taken into account an extra idle time of 
magnitude ^Tidio = ti" after every atom crosses the cavity. This 
might be necessary to ensure the right spacing between them. 
As a result, one 3-qubit graph state is yielded every 30.2 /is 
and 62.6 /is on average, respectively. We have also studied 
the resilience of our protocol versus a hypothetical fluctuation 
in the interaction time. This also models experimental imper- 
fections in selecting atomic velocities and controlling the de- 
tunings. Figs. [3](c) and (d) show the decrease in fidelity and 
probability of heralding, respectively, assuming that there is 
up to a 10% uncertainty (uniformly sampled) in the selection 
of the interaction time. The numerical results indicate that our 
scheme is robust against this type of imperfection. 

As a final remark, the most crucial aspect to judge if the 
generation of the graph state has been successful after all is 
to determine whether the state j^a) is really entangled or not. 
In spite of the lack of general criteria for multipartite entan- 
glement, we can still make use of the entanglement witness 
derived in Ref. [2^^9J- In our case, this particular type of 
witness for graph states allows to relate the existence of en- 
tanglement in the state with its fidelity. More precisely, the 
witness says that any state close to a graph state with fidelity 
larger than 1/2 should be entangled. 
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